European Journal of Histochemistry 2013; volume 57:e24 



Correspondence: Federico Boschi, Department of 
Neurological, Neuropsychological, Morphological 
and Motor Sciences, University of Verona, strada 
Le Grazie 8, 37134 Verona, Italy. 
E-mail: federico.boschi@univr.it 

Key words: Lipid droplet, 3T3-L1, adipocyte, fat, 
triglyceride accumulation, integrated optical den- 
sity. 

Acknowledgments: the authors would like to 
thanks Dr. Lucia Calciano for her useful discus- 
sions in the statistical field. 

Contributions: VR, data acquisition, analysis and 
interpretation; FB, concept, analysis and inter- 
pretation of data, draft of the article; MP, EZ, data 
acquisition; AS, MZ, revision of article for impor- 
tant intellectual content. 

Conflict of interests: the authors declare no con- 
flict of interests. 

Received for publication: 27 May 2013. 
Accepted for publication: 27 June 2013. 

This work is licensed under a Creative Commons 
Attribution Noncommercial 3.0 License (CC BY- 
NC 3.0). 

©Copyright V. Rizzatti et ai, 2013 

Licensee PAGEPress, Italy 

European Journal of Histochemistry 2013; 57:e24 

doi:10.4081/ejh.2013.e24 



Lipid droplets characterization 
in adipocyte differentiated 
3T3-L1 cells: size and optical 
density distribution 

V. Rizzatti, 1 F. Boschi, 2 M. Pedrotti, 1 
E. Zoico, 1 A. Sbarbati, 2 M. Zamboni 1 

'Department of Medicine, Geriatric 
Section, University of Verona 
department of Neurological, 
Neuropsychological, Morphological and 
Motor Sciences, University of Verona, 
Italy 



Abstract 

The 3T3-L1 cell line, derived from 3T3 cells, 
is widely used in biological research on adipose 
tissue. 3T3-L1 cells have a fibroblast-like mor- 
phology, but, under appropriate conditions, they 
differentiate into an adipocyte-like phenotype. 
During the differentiation process, 3T3-L1 cells 
increase the synthesis of triglycerides and 
acquire the behavior of adipose cells. In partic- 
ular, triglycerides accumulate in lipid droplets 
(LDs) embedded in the cytoplasm. The number 
and the size distribution of the LDs is often cor- 
related with obesity and many other patholo- 
gies linked with fat accumulation. The integrat- 
ed optical density (IOD) of the LDs is related 
with the amount of triglycerides in the droplets. 
The aim of this study is the attempt to charac- 
terize the size distribution and the IOD of the 
LDs in 3T3-L1 differentiated cells. The cells 
were differentiated into adipocytes for 5 days 
with a standard procedure, stained with Oil Red 
O and observed with an optical microscope. The 
diameter, area, optical density of the LDs were 
measured. We found an asymmetry of the ker- 
nel density distribution of the maximum Feret's 
diameter of the LDs with a tail due to very large 
LDs. More information regarding the birth of 
the LDs could help in finding the best mathe- 
matical model in order to analyze fat accumula- 
tion in adipocytes. 



Introduction 

The role of fat accumulation has emerged as 
a fundamental target of medical studies 
because it is responsible of obesity, a wide- 
spread disease in affluent societies. Obesity is 
a major cause of preventable death in the 
world. 13 In particular, excessive triglycerides 
accumulation in humans can lead to cardiovas- 
cular dysfunction, diabetes mellitus type 2, 
hypertension, hypercholesterolemia, metabolic 



syndromes, hypertriglyceridemia, disorders of 
osteo-articular system, stroke, sleep apnea syn- 
drome and some types of cancer. 45 In recent 
years, in order to understand the molecular 
regulation of human obesity, immortalized 
murine cell lines such as 3T3-L1 cells and 
F442A were studied. 6 3T3-L1 is a cell line 
derived from 3T3 cells. 3T3-L1 cells have a 
fibroblast-like morphology, but, under appropri- 
ate conditions, they differentiate into an 
adipocyte-like phenotype The differentiation 
process of 3T3-L1 cells in adipocyte-like cells is 
based on the increased synthesis and accumu- 
lation of triglycerides, for this reason the 3T3- 
Ll cells are extensively used to study adipogen- 
esis. 7 Triglycerides accumulate in lipid droplets 
(LDs) in the cytoplasm. LDs are intracellular 
organelles that regulate the storage and hydrol- 
ysis of neutral lipids, including triglycerides 
and cholesterol esters. 8 

The Oil red O staining was used in many dif- 
ferent work to quantify the adipose conversion 
and triglycerides 9 and loss of lipids in decellu- 
larized tissue. 10 The integrated optical density 
(IOD) of the LDs was used to quantity the lipid 
accumulation in muscle fibers which was found 
significantly increased in obesity and reduced 
with weight loss. 11 Moreover, Ramirez-Zacarias 
and co-workers found that, staining the 
adipocyte with Oil red O, the IOD increases lin- 
early with the triglycerides content in 3T3 
adipocytes. 9 LDs are dynamic and heteroge- 
neous in size, location, and protein content. In 
literature it has been described the cellular role, 
processes of biogenesis, mobilization of fat and 
role of proteins associated with LDs 12 but no 
one has comprehensively described the impor- 
tant aspect of their heterogeneity. The aim of 
this work was to quantify the heterogeneity in 
size and optical density of intracellular LDs in 
3T3-L1 cells. This represents the first step in 
order to analyze fat accumulation in cells and 
the changes in pathological disorders. 



Materials and Methods 

Cell line and cell culture 

The 3T3-L1 preadipocytes (ECACC, Sigma- 
Aldrich, St. Louis, MO, USA) were cultured in 
250 mL polystyrene sterilized flasks equipped 
with 0.2 urn vented plug seal cap by BD (Becton 
Dickinson, Franklin Lakes, NJ, USA), with 
DMEM/GlutaMAX, supplemented with 10% of 
FBS, and 1% of Antibiotic Antimycotic Solution, 
to reach a concentration of 2.5-3.0xl0 6 cells 
(confluence 85-90%). Then, 3T3-L1 cells were 
detached from flask with trypsin-EDTA and fur- 
ther eluted to be seeded in 4-wells BD Bio 
Coat™ CultureSlides. Two days after conflu- 
ence in wells, cells were washed three times 



with PBS and induced to differentiation with 
DMEM/F12 containing 10% FBS, 1% Antibiotic 
Antimycotic Solution, 0.2 mM IBMX, 10 \M 
rosiglitazone, 1 uM dexamethasone, 10 ^ig/mL 
insulin for three days. After 72 h /?a#-induction, 
cells underwent PBS washing (three times), 
and medium was replaced with an adipocyte 
maintaining medium (AMM), composed of 
DMEM/F12 enriched with 10% FBS, 1% 
Antibiotic Antimycotic Solution, 10 pig/mL 
insulin, in which cells were cultured for 2 days. 
3T3-L1 cells were incubated at 37°C, 5% C0 2 in 
each culture step. 

Oil Red O staining 

After washing with PBS 0.1M pH 7.4, cell cul- 
tures were fixed for 20 min with 4% formalin in 
PBS 0.05M. Cells were washed with sterile dou- 
ble distilled water and subsequently with 60% 
isopropanol for 2 min and stained with a filtered 
0.35% Oil Red O solution in 60% isopropanol for 
10 min at room temperature. Then, cells were 
washed with sterile double distilled water and 
stained with Mayer's Hematoxylin Bio-Optica 
ready to use solution for 1 min at room temper- 
ature and washed again with sterile double dis- 
tilled water. Slides were treated with Dako fara- 
mount aqueous mounting medium ready to use 
and then was applied to the coverslip. 
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Image and statistical analysis 

Cells were observed in an Olympus BX51 
photomicroscope equipped with a KY-F58 CCD 
camera (JVC). The images were analyze and 
stored using the software Image-ProPlus on a 
personal computer. Ten cells were imaged and 
1935 LDs were detected. For each LD area, 
perimeter, Maximum Feret's Diameter (MFD) 
and IOD were measured. The Feret's diameter 
(FD) is a measure of an object size along a 
specified direction, defined as the distance 
between the two parallel planes (or lines in 
case of planar images) restricting the object 
perpendicular to that direction. FD is generally 
used in optical microcopy to measure the size 
of irregularly shaped particles. 13 Here we 
measured the MFD, the longest distance 
between any two points of the object, also 
known as maximum caliper diameter. The LDs' 
properties were measured using a semi-auto- 
matic procedure based on the automatic rou- 
tines implemented in ImageJ with a strict con- 
trol of an operator. The procedure and the 
parameters used to guide the routines were 
previously carefully checked and the results 
were compared with the ones obtained com- 
pletely manually from an operator on a subset 
of LDs (-200). With the help of contrast varia- 
tion and color changing the operator was able 
to separate almost all the clusters of LDs in 
single LDs. The non resolved clusters were 
removed from the analysis, they were very few 
with respect the total amount of LDs (less than 
3.3%). The statistical analysis was conducted 
with both Matlab R2007a (Mathworks) and 
EasyFit 5.5 (Mathwave). 



Results 

Size distribution 

The images show that LDs are well stained 
and the image contrast is good (Figure la). 
This result is useful for the semi-automatic 
analysis procedure. The form of the 1935 LDs 
detected results almost perfect circular. 
Plotting area versus perimeter (Figure lb) this 
visual impression is confirmed by the distribu- 
tion of the experimental data in comparison 
with the perfect circularity represented by the 
red line. 

The perfect circularity is verified when Area 
= 4jtx [Perimeter] 2 . So the experimental data 
were fitted with a quadratic curve and the fit 
(dotted line) is very good (R2=0.999). 

Figure lc shows the histogram of the MFD 
of the LDs. LDs have MFD in the range 0.3-3.6 
^im. The histogram is clearly asymmetric and 
shows a long tail in the positive axis reflecting 
the presence of LDs with sizes very higher 
with respect the mean of the population. To 



avoid the dependence of the histogram repre- 
sentation to the choice of the sub-intervals or 
bins which cover the range of the data, the ker- 
nel density of the MFD distribution was calcu- 
lated and it is visible in Figure Id. Most of the 
LDs has a MFD around 0.5 urn. The more 
detailed representation of the data suggests a 
local secondary peak around 0.8 ^im. 

In order to analyze the distribution of the 
MFD of the LDs we tested the goodness of dif- 
ferent distributions. In the widespread of 
asymmetric distributions we found that the 
experimental data can be fitted with a subset 
of distribution, most of them described with 
more than three parameters. In order to sim- 
plify the model we focused on distributions 
with no more than three parameters; in partic- 
ular we found that the data can be fitted with 
the Generalized Extreme Value (GEV), lognor- 
mal and gamma functions. The GEV distribu- 
tion is a family of continuous probability distri- 
butions developed within extreme value theory 
to combine the Gumbel, Frechet and Weibull 
families also known as type I, II and III extreme 
value distributions (defined by three parame- 
ters: local parameter u, scale parameter o, 
and shape parameter k). The lognormal distri- 
bution is a continuous probability distribution 
of a random variable whose logarithm is nor- 
mally distributed (defined by two parameters: 



a 




the shape parameter o and the log-scale 
parameter u). The gamma distribution is a 
two-parameter family of continuous probabili- 
ty distributions with different parameteriza- 
tions: we used here the shape parameter a and 
the rate parameter |3. We obtained the subse- 
quent fitting values: k=0.15731, a=0.27612 
and u=0.64964 for GEV, a=0.46083 and 
p=— 0.26163 for lognormal and a=3.9329 and 
|3=0.21852 for gamma distribution. The good- 
ness of the fitting was evaluated with tree dif- 
ferent methods: GEV distribution obtained the 
fit with statistical test Kolmogorov Smirnov (K- 
S)=0.16886 and rank 34, Anderson Darling (A- 
D)=33.3222 and rank 20, x 2 =1809.3 and rank 
43. Lognormal distribution obtained: statistical 
test K-S=0.16661 and rank 31, A-D=32.865 and 
rank 18, x 2 =1714.5 and rank 37. Gamma distri- 
bution obtained: statistical test K-S=0.13633 
and rank 2, A-D=35.777 and rank 25, x 2 
=382.59 and rank 1. Despite all the three 
distributions can fit the data the gamma distri- 
bution obtained the best fit (lower rank 
values). The experimental data and the three 
distributions are visible in the same plot in 
Figure 2. Despite GEV distribution fits the 
peak of the experimental data distribution bet- 
ter than the others distributions, the gamma 
distribution is the best fitting curve conside- 
ring the entire set of data. 





Figure 1. Oil Red O stained lipid droplets are well visible in the peri-nuclear region of a 
3T3-L1 cell differentiated into adipocyte (a). The circularity of the lipid droplets is evi- 
dent in the plot area vs perimeter (b); the red line represents the perfect circularity and 
the dotted line the best fit. Histogram (c) and kernel distribution (d) of the maximum 
Feret's diameter of the lipid droplets. 
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Optical density distribution 

Figure 3a shows the histogram of the IOD of 
the LDs. Most of the LDs has low IOD, instead 
a small number presents a higher IOD as 
reflected in the asymmetric distribution. The 
kernel distribution is visible in Figure 2b. The 
values span from 2 to 410. In order to fit the 
experimental data distribution we tried with 
the asymmetric curves described above: GEV, 
lognormal and gamma distributions. We 
obtained the subsequent fitting values: 
k=0.27225, o=15.93 and u=25.809 for GEV; 
a=0.6653 and u=3.4766 for lognormal; 
a=1.5249 and (3=26.755 for gamma distribu- 
tion. Also for IOD the goodness of the fitting 
was evaluated with tree different methods. 
GEV distribution obtained the fit with statisti- 
cal test K-S=0.03519 and rank 18, A-D=2.7728 
and rank 15, x 2 =56.292 and rank 18. 
Lognormal distribution obtained: statistical 
test K-S =0.03207 and rank 13, A-D=2.4218 
and rank 11, x 2 =48.446 and rank 15. Gamma 
distribution obtained: statistical test K- 
S=0.12422 and rank 30, A-D=42.76 and rank 
27, x 2 =299.27 and rank 29. Despite all the 
three distributions can fit the data only the log- 
normal distribution obtained the best fit. The 
three distributions used for the fit and the 
experimental data are visible in Figure 3c. 



Discussion 

In recent years ultrastructural analysis 
revealed the morphological polymorphism of 



the adipose tissue with different lipid 
depots. 1415 This results encourages the study of 
lipid storage in adipocytes. This work demon- 
strates that 3T3-L1 cells represent an excellent 
model to study the process of lipid accumula- 
tion in cytoplasm. The shape of the detected 
LDs resulted almost perfectly circular suggest- 
ing a homogeneous distribution of the trans- 
port through the surface of the LDs. 



Figure 3. Histogram (a) 
and kernel distribution 

(b) of integrated opti- 
cal density of the lipid 
droplets (LDs). Com- 
parison between kernel 
density of the integrat- 
ed optical density data 
and three fitting curves: 
Generalized Extreme 
Value, lognormal and 
gamma distributions 

(c) . The lognormal dis- 
tribution represents the 
best fit of the data. 




Figure 2. Kernel density of the maximum Ferets size of the lipid 
droplets and three fitting curves: Generalized Extreme Value, log- 
normal and gamma distributions. Despite Generalized Extreme 
Value distribution well fits the peak of the distribution, gamma 
distribution represents the best fit with respect the entire dataset. 



The planar images of the differentiated 3T3- 
Ll cells showed some overlaps between the 
LDs. Many of them can be resolved using dif- 
ferent thresholds and varying the contrast 
parameters in the representation of the 
images. The non-resolved clusters of LDs are 
very few and can be neglected when using a 
great number of LDs. In the planar images it is 
possible that small LDs located below the 
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largest LDs along the z-axis (respect the cells 
layered in the xy plane) were not observed. A 
new interesting approach could derive from the 
use of the confocal microscopy to analyze 3D 
images to extract a more detailed size distribu- 
tion of the LDs. We described here the MFD and 
IOD distributions of the LDs in 5-days differen- 
tiated 3T3-L1 cells. The continuous distribution 
of the IOD shows the absence of eventual sub- 
populations of LDs characterized by different 
amounts of lipids. Ramirez-Zacarias and co- 
workers found that Oil red 0 specifically stains 
triglycerides and cholesteryl oletate ,but no 
other lipids and showed a linear relationship 
between the degree of staining and the lipid 
contents in fat cells. 9 In case of 3T3 adipocytes 
lipid accumulation is mainly due to increase of 
triglyceride synthesis, but cholesterol or its 
derivatives did not seem to accumulate during 
adipocyte differentiation. 16 Thus, it is possible 
to argue that the IOD curve presented here 
reflects also the continuous distribution in the 
triglycerides content in the LDs. The IOD data 
represent a template in order to test the effec- 
tiveness of different treatments. In particular 
variations in the IOD data distribution could be 
expected when the cells increase the synthesis 
of the triglycerides or the concentration of the 
lipids in the LDs. Comparing the IOD distribu- 
tion before and after treatment it could be pos- 
sible to acquire information on the lipids con- 
centration in the new born LDs. 

Many different pathologies involved in fat 
accumulation are related to changing in LDs 
number and size. Some studies analyzed the 
size of the LDs and the changing in the size 
after different treatments but they focused on 
the mean value and the standard deviation 17 or 
on a subset of LDs randomly selected. 18 In this 
study we measured all the detectable LDs in 
the analyzed cells and we propose a mathemat- 
ical description of the entire size distribution 
of the LDS in order to create a standard popu- 
lation for future comparisons in different 
(healthy and pathological) conditions. The 
description of the entire size distribution can 
overcome the limitation due to the comparison 
between only the mean values and, for exam- 
ple, it can show differences in distributions 
with the same mean values. In general the 
curve representing a distribution contains 
many information not well represented by a 
single value. In this study we used curves char- 
acterized by two or three parameters which are 
related to the mean, the variance, the skew- 
ness and other properties of the distribution 
itself and a description with more than a single 
parameter could highlight changes, for 
instance, in the tail of the distributions. 

We found that the best distribution to fit the 
distribution of the size of the LDs is a gamma 
curve. The local maximum in the MFD distribu- 
tion around 0.8 urn represents a problem for 



the fitting. We hypothesized that it could derive 
from the quantization of the measurements 
done on the computer screen. Using different 
magnifications of the images at the monitor on 
a sub-population of LDs (-200) we did not 
observe the peak. Otherwise, the local maxi- 
mum could be due to non homogeneity of the 
analyzed cells; in particular differences in LDs' 
size could reflect different ages of the cells 
themselves because it was observed that the 
size of the LDs can change along time. 19 
Another hypothesis is that the entire distribu- 
tion could be considered the sum of two differ- 
ent distributions; in this case the increase in 
size of the LDs could occur in two or more 
steps: synthesis of triglycerides and formation 
of the small LDs (first step) and than fusion of 
the LDs (second step). The asymmetry of the 
kernel density distribution of the MFD of the 
LDs could be related with the increasing of the 
sizes of the LDs along time due to the fusion of 
the LDs themselves. Fusion of LDs was effec- 
tively observed by different authors with confo- 
cal microscopy. 19 20 We found the best distribu- 
tion for MFD and IOD of the LDs, but at the 
present time more than a single model could be 
valid to describe the LDs in 3T3-L1 5-days dif- 
ferentiated adipocytes. A mathematical de- 
scription or more information regarding the 
birth of the LDs could help in finding the best 
mathematical model. In any case, we would like 
to propose a mathematical approach to quanti- 
fy the modifications in order to evaluate many 
different pathologies or the growth of the LDs 
along time. 
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